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We have examined the physics and the experimental feasibility of studying various kaon 
decay processes in which the polarization of a muon in the final state is measured. 
Valuable information on CP violation, the quark mixing (CKM) matrix, and new physics 
can be obtained from such measurements. We have considered muon polarization in 
— > n'^ tJ.~ and K — ► ■jTfi'^ fi~ decays. Although the effects are small, or difficult to 
measure because of the small branching ratios involved, these studies could provide clean 
measurements of the CKM parameters. The experimental difficulty appears comparable 
to the observation of A" — > ttui/. New sources of physics, involving non-standard CP 
violation, could produce effects observable in these measurements. Limits from new 
results on the neutron and electron electric dipole moment, and ^ in neutral kaon 
decays, do not eliminate certain models that could contribute to the signal. A detailed 
examination of muon polarization out of the decay plane in /x+tt"!/ and — > 

fi'^i/^ decays also appears to be of interest. With current kaon beams and detector 
techniques, it is possible to measure the T-violating polarization for A'+ — > /z+tt";/ with 
uncertainties approaching ~ 10~^. This level of sensitivity would provide an interesting 
probe of new physics. 

1. Introduction 

We have examined the possibihty of measuring muon polarization asymmetries that 
are sensitive to P, T or CP symmetries; these are tabulated in Table |^. Observation 
of the possible effects requires high fluxes (> 10^^ K decays per year) of kaons, now 
available at several accelerator facilities, notably at Brookhaven National Labora- 
tory AGS and at the KEK-PS. In the near future, the Fermilab main injector, as well 
as the Japanese Hadron Factory, could deliver much higher intensities of separated 
as well as unseparated kaon beams. The 4> meson factory, DAPHNE, is also a new 
facility with an intense pure-kaon flux. Thus far, CP violation has been observed 
conclusively only in the neutral kaon system. Although a theoretical description 
of CP-violation in the neutral kaon system is available in the single complex phase 
of the standard-model Cabibbo-Kobayashi-Maskawa (CKM) matrix, part of, or the 
entire phase, can have origin in deeper causes that have so far eluded experimental 
scrutiny. During this past decade, experiments at FNAL and at CERN, focusing 
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Table 1: Decay modes and polarization asymmetries (or correlations) of interest in 
K decays. The symbols s and p refer to the spin and momentum vectors in the 
decays. 



on the measurement of the direct tttt transition, or — , have reported ever- 

improved results, the latest being Re(e7e) = (28.0 ± 3.0 ±'2.8) x lO"'^ (FNAlS) 
and Re(e7e) = (14.0 ± 4.3) x lO""* (CERnI). These provide conclusive evidence 
for the presence of direct CP violation in K'L decays, and there is great theoretical 
effort in progress to interpret these numbers.0~B We will not attempt here to review 
e'/e; but it is clear that it is not yet certain if the standard model description of 
CP violation can accommodate these results. In addition to e'/e, rare kaon decays 
are also of interest for understanding the CKM matrix. For a recent review of rare 
kaon decay processes see Ref. |^. 

Over the next decade, ambitious efforts towards gaining a better understanding 
CP-violation and the CKM matrix elements will be pursued at B-factories. The 
importance of these efforts is undeniable, but it is also worthwhile to investigate 
the possibility that some or all of the CP-violation arises from effects outside of the 
minimal standard model, particularly outside of the current CKM matrix. 

It should be recalled that CP-violation is required to generate the observed 
baryon asymmetry in the universe, and it is now accepted that the CP-violation em- 
bodied in the CKM matrix does not have sufficient strength to serve this purposeJl 
It is important to examine if sources of physics beyond the standard model that 
could generate the baryon asymmetry can also generate CP or T violating muon 
polarizations in the kaon decay modes given in Table 0. 

We will now briefly consider the measurement of muon polarization. This has 
been discussed before in the literature.Blla Muon polarization in the experiments 
under consideration would be observed by stopping the muons in some appropriate 
material, and measuring the direction of the decay electron (positron in the case of 
The muon decay spectrum is given by: 

(3-2z)±|P|cos6'(2z-l)) (1) 



dzdQ 2n 

where the positive sign in the brackets is for /i+ decay and the negative sign is for 
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Mode Branch. Standard Final Non-SM Rcf. 

Fraction Model State Int. value 

(1) K+ -> TT^'n+iy 0.032 0.0 ^ 10"*^ < 10'^ 

(2) K+ jjL+v-j 5.5 X 10--^ 0.0 ^ 10'^ < 10"^ 

(3) Kl n+fj- 7.2 X 10-" ^ 0.002 0.0 < 10"^ 
~lt) K+ TT+n+n- 7.6 X 10-» - 10-^ ~ = 

(5) 0.0 - 10-3 ~ 10-3 

(6) ^0.06 ~0.0 ^0.1 

~(7) Kl TT^'fi+fi- - 5 X 10-1^ = - - 

(8) ~ 0.5 - - 

Table 2: The decay modes and asymmetries of interest; the row numbers corre- 
spond to those in Table ^ The other columns are: the known standard-model 
branching ratio (note that the Kl — > t:^ ^'^ ^~ branching ratio is not yet measured, 
the present 90% C.L. limit is < 3.8 x 10~^^), the estimated standard-model value 
of the asymmetry, the value due to final-state interactions, the maximum possible 
value allowed by non-standard physics, and the theoretical reference. Some of the 
results have been adjusted to account for more recent values of the top quark mass 
(174 GeV/c^). In the case of K^ /^^/^^ ^nd K^ jjt^ ji^ , the theoretical 

estimates for maximum possible non-standard contributions to T violation do not 
agree; we chose the mean value of the different estimates. The "-" means that there 
is no reliable prediction. 



fj,~ decay, z — 2Ee/mf^, and 9 is the angle between the polarization P and the 
direction of the positron or electron. We will now restrict ourselves to only /i+ 
decays. The magnitude of the polarization P, multiplied by a dilution factor D(zo), 
is the asymmetry A in the number of decays that produce positrons forwards versus 
backwards with respect to a plane normal to the polarization vector. 

where A'^i and N2 are the number of forward and backward decays, respectively. 
D{zq) is the asymmetry dilution factor that depends on the lower-energy cutoff 
zq = 2E{)/m^, where Eq is the minimum observable energy of the electron from 

jjL decays. The uncertainty on the asymmetry is given by 5A = \J 1 where 
TV = Ni -\- N2. When the asymmetry is small, the error on the polarization is 
given by SP ~ practical devices, there is usually a lower cutoff on 

the positron energy. Integrated over the entire spectrum, the asymmetry is |^'|/3. 
Clearly, both the total number of detected decays N, and the dilution factor D, 
depend on the low-energy cutoff. Best performance is reached when the lower 
cutoff is at about zq 0.75, which keeps approximately 1/2 of the spectrum and 
corresponds to a dilution factor of 0.5. Often, the muons are precessed by a small 
magnetic field (at a rate of 42.5 kHz per Gauss) perpendicular to the direction of 
the spin. This makes it possible to measure two components of the polarization, as 
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well as eliminate systematic differences in detector efficiencies, at only a cost of \/2 
loss in the dilution factor. 

Other considerations that determine the asymmetry are depolarization of the 
stopped /x"*", and confusion from the multiple scattering of the decay positron. It 
is well known that multiple Coulomb scattering of the muon, as it comes to rest, 
does not cause depolarization; however, after it slows down to atomic velocities, 
there are many processes that can lead_to depolarization. The amount of depolar- 
ization is strongly material dependent Carbon (in graphite form) and aluminum 
are considered good materials for muon polarimeters because they preserve the po- 
larization. Detectors such as wire chambers or scintillators must be placed next to 
blocks of graphite or aluminum to detect the decay positrons. No practical active 
materials such as plastic scintillator or scintillating crystals have been found that 
would preserve polarization. Dilution of the asymmetry from multiple scattering 
of the positron depends on the geometry and the energy cutoff. Ultimately, the 
dilution factor must be measured in each experiment to determine the sensitivity 
of the result. 

2. K+ n^n+iy 

The transverse, or out-of-plane, muon polarization in this decay has been analyzed 
previouslyJljila The decay amplitude, Ai, can be written as follows: 

M = ^ sine J+{q') {{PK+P.f + ^{q^){pK ~ P.^)) {u^,lx{l ~ (3) 

where Gp is the Fermi constant, sin^c is the Cabibbo angle, and q, pK, Ptt are 
the momentum transfer, the kaon and the pion 4-momenta, respectively. The out- 
of-plane polarization (Pj") of the muon is non-zero when the form factor ^ has an 
imaginary component. This polarization is expected to vanish to first order in the 
standard model, because of the absence of the CKM phase in the decay amplitude. 
Irreducible backgrounds, such as from final-state intera,ction s ( FSI) in this decay 
are expected to be small (~ 4 x 10~^), and can be ignored.c§E3 It has been shown 
that any model involving only effective V or A interactions cannot produce this 
type of polarization.EZI The existence of a non-zero value of this polarization will 
therefore provide a definite signature of physics beyond the standard model. Pj' is 
a function of the two Dalitz variables that define the 3-body — > /i+tt^i/ decay. 
In most experimental situations, one averages over some portion of the Dalitz plot. 
Therefore, the average has two components: (1) a kinematic factor {fo) that 
describes the average over the Dalitz plot, including experimental acceptance, and 
(2) the imaginary part of an amplitude or combination of amplitudes (Im^). We 
will often use the same P^ to denote the average over an ensemble of events: 
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/d is estimated by setting ^ = 0, and ignoring terms of O (—-2-\ 

/ m„ sin^M^ \ 



JD 



ruK +COS V, 

I -^M I 



(4) 



where is the muon energy and 0^^ is the angle between the muon and the 
neutrino in the kaon rest frame. The dependence over the Dahtz plot is averaged 
with the experimental acceptance. If we choose a typical point on the Dalitz plot, 
e.g., K, nOMeV and 6^^ = 90°, then the value of fn calculated at that point 
is fo ~ 0.17. The value of Im^ is model dependent; for example, in the case of a 
non-standard effective scalar interaction, Im^ is proportional to the imaginary part 
of the scalar coupling strength. 

In particular, multi-Higgs and leptoquark models can produce non-zero out-of- 
plane polarization. In multi-Higgs models, a charged Higgs particle mediates a 
scalar interaction that interferes with the standard model decay amplitude; in such 
models, the polarization can be as large as 10"'^, without conflicting with other 
experimental constraints, for example, the measurements of the neutron electric 
dipole moment and the branching fractions for B Xtv, and b — > 57. The 
indirect limits on Pj" from other measurements ha v e been examined in the context 
of the minimal 3 Higgs Doublet Model (3HDM) .EllO'ca The transverse polarization 
for 3HDM is given by 

2 2 I 2 I 2 tr\ 

V = v^^v^+v^ (5) 

where s^, S3, and C3 are unknown parameters from the unitary 3x3 mixing matrix 
of the 3 Higgs doublets. It is assumed that the ratios of the 3 vacuum expectation 
values are the same as the ratios of the 3 generation masses; v\ : V2 vz Wb : 
mt : rxir- We have updated the calculation of constraints on the 3HDM model to 
include new experimental results and collected them in Table |[ The constraints are 
calculated in terms of two assumptions on the mass of the charged Higgs: m;i ~ mw 
or m/i ~ 2mw- Both of these assumptions are above the current lower limit of m/j 
78GeV on the mass of a charged Higgs boson from LEP.cJ The best constraints are 
from the branching ratio measurements B{b — s- 57) and B{b — s- Xtv). These are 
unlikely to improve in the near future because the errors have large theoretical 
components. Furthermore, these constraints have the requirement that the real 
part of the 3HDM amplitude cancels with the standard model amplitude. Without 
this assumption, the constraints would be quite weak. Nevertheless, as can be seen 
in the table, even these optimistic limits allow any value of below the current 
direct bound, which is described below. T-violation could occur in the minimal 
supersymmetric models through an effective scalar interaction involving the charged 
Higgs particle, nevertheless, the effect isixuisidered too small to observe, except in 
models that contain R-parity violation.! 
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Measurement 


Value 


< 

rrih ~ mw rrih ~ 2mw 


dn 
de 

it 

niK 
B{1 
B{1 


34] 


< 7.5 X lO-^^e - cm (95%C.L.) 

< 1.6 X 10^26g _ (95%C.L.) 
(1.9 ±0.24) X 10-3 

(0.5301 ±0.0014) X 10i°?is-i 
(2.54 ±0.56) X lO-'^ 
(4.08 ±0.98) X 10-2 


0.039 0.039 

0.09 0.02 
0.009 0.009 



Table 3: Ccpstraints on P^ of ^^-k^v decay for 3HDM from other 

measurementstM3 The "-" means that there is no significant constraint. AU mea- 
surements have been scaled to correspond to the same 95% confidence level to get 
consistent constraints. 



The best lirmt on this process was obtained recently by an experiment at the 
KEK-PS, E246.ii They measured Pj = -0.0042 ± 0.0049 ± 0.0009, which corre- 
sponds to a value of the T violating parameter: Ini^ = —0.013 ± 0.016 ± 0.003, 
where the errors are statistical and systematic, respectively. This measurement 
was performed using approximately 3.9 ^f-XO^ events. Additional data may in- 
crease the total sample by a factor of two.cHi Experimental limits were obtained al- 
most 20 years ago with both neutral and charged kaons at the BNL-AGS.Ol3 The 
experiment with decays made a measurement of the transverse polarization, 
= 0.0031 ± 0.0053. The combination of the neutral and charged experiments 
could be interpreted as a limit on Im^ = —0.01 ± 0.019. 

Measurements from the 1980s were based on 1.2 x 10'' and 2.1 x 10*" K'^ 
decays to fj,~^TTi>, and were limited by low analyzing power and by backgrounds. 
The KEK-E246 measurement relied on the new technique of using a stopped if + 
beam, and measuring the muon decay direction in an aluminum absorber, without 
spin precession. The BNL measurements used in-flight kaon decays, stopped the 
muons in an aluminum absorber, and measured the muon spin by its precession 
in a weak magnetic field. Both techniques relied on the cylindrical symmetry of 
the apparatus to suppress systematic errors due to non-uniform efficiencies. The 
systematic error in KEK-E246 was reduced further by collecting events with for- 
ward and backward going pions. The BNL experiment alternated the direction 
of the precessing magnetic field to cancel the detection efficiency to second order. 
The remaining systematic errors in both techniques are quite different: the largest 
systematic error in the stopped kaon experiment comes from the knowledge of the 
fringe magnetic field in the muon stopping region. The largest systematic error for 
the in-flight experiment is thought to be due to the mechanical misalignment of the 
azimuthal segments of the muon absorber with respect to the precessing axis of the 
magnetic field. 

A new experiment should reach much higher statistics and have high analyzing 
power for the muon polarization. KEK-246 has demonstrated that a well designed 
muon stopper with low background can have an analyzing power of 0.197 ± 0.005. 
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It is possible Id extend the sensitivity of the KEK-246 experiment with a more 
intense beam,E3 such as the low energy separated beam used by BNL experiment 
E787,lJ which has been able to obtain kaon stopping rates of approximately 1 million 
per sec. A stopped kaon experiment, with a field for precession has the potential 
advantages of having better systematic control using cancelations frorn both the 
forward/backward 7r° direction, as well as by alternating field direction.c3 

A new experiment has been designed at the BNL- AGS to perform this measure- 
ment, with an error on the polarization approaching 10"** .El The design is based on 
the 1980 experiment, however, it uses a 2 GeV/c separated charged kaon beam to 
reduce the background counting rate. Other improvements will involve higher ac- 
ceptance and analyzing power, with a larger apparatus and a more finely segmented 
polarimeter made of graphite. The experiment will collect approximately 550 events 
per AGS pulse (3.6 seconds). Thus the statistical accuracy of the polarization mea- 
surement in a 2000 hr (2 x 10^ pulses) run could be 5Pt ~ 1-3 x 10~^. Through a 
precise alignment of the polarimeter and through measurements of null asymmetries 
from ~f 11^1/ decays the experimenters expect to control the systematic error 
in the apparatus to below the statistical uncertainty. Other studies carried out for 
experiments at the Japanese Hadron Factory and the (j) factory a|LllAPIINE also 
show that statistical sensitivities of the order ~ 10"^ are possible.cZrEa 

Finally, we note that it could be easier to measure the transverse polarization 
in the case of 'k~ fx^v^ decays, because it is easier to detect a charged pion. 

However, the transverse polarization in this case is contaminated by final state 
interactions due to the presence of a charged pion. This effect, averaged over the 
phase space, has recently been estimated to be = 2.4 ± 0.1 x 10~'^.CJ 

3. K+ ^+1/7 

The branching ratio for decay is (5.5 ± 0.28) x 10^^.0 The decay 

is dominated by inner bremsstrahlung. Recently, the structure-dependent part of 
the branching ratio (mostly for positive-helicity photons) has been measured to be 
(1.33 ± 0.12 ± 0.18) X IO^^.IhI The structure-dependent form factors have also been 
measured recently in decays of if+ e+i/e+e^ and if+ /i+i/e+e~, and are 
sensitive to the form factor for negative-helicity photons.0 

In fj.'^v"/ decay, the transverse muon polarization, which is T-violating, 

can arise from interference between inner bremsstrahlung and the structure-dependent 
part of the decay. It is sensitive to new pseudo-scalar, vector, and axial- vector intet 
actions, in contrast to — > /z+tt^^ decay's sensitivity to new scalar interactions.Ej 
Recent studies have shown that in certain extensions of the standard model, includ- 
ing SUSY, there coul d be a T-violating muon polarization in /i+i^7 decays 
as large as 10~^.E3CJ Others argue that it is unnatural to generate a contribution 
to transverse polarization in K'^ ^~^vj decay larger than 10~^ in SUSY, unless 
R-parity is brokenO Finally, final-state interactions can induce a non-zero trans- 
verse polarization. This haabeen calculated recently to be of the order of 10"'', and 
varies on the Dalitz plot.oEj It can be as large as 5 x 10~^ at the high end of the 
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muon energy spectrum, for ^ Mk/4:. 

Following the calculations in Ref. the K'^ ^'^vj decay is described by 
"inner bremsstrahlung" (IB) and "structure-dependent" (SD) contributions. The 
parameters are the kaon decay constant fx, axial vector form factor Fj^, and a 
vector form factor Fy ■ Physics beyond the SM is introduced in additional terms in 
the Lagrangian 

Gf _ 

C = ---^si7i6'cS7"(l-75)-uz^7a(l-75)^ + Gssuj^(l+75)^ + GpS75Ui/(l + 75)A* 

+ GyS7"-ui?7a(l - 75)/^ + GAS7"75Ui^7a(l - 75)m + ^-c-, (6) 

where Gf is the Fermi constant, 9c is the Cabibbo mixing angle, and Gs, Gp, Gy, 
and Ga are parameters of the new interactions due to scalar, pseudo-scalar, vector, 
and axial vector exchange, respectively. 

The new interactions will modify the decay constants and form factors in the 
following way: 

fK = f%{l + Ap + AA), 

Fa = F^(1 + Aa), 

Fv = FO(I-Av), (7) 

and 

A ^ f_Gprr^^ r r \ 

^iP,A,v) - 7; 7 ] ^ ,G^,Gy . (6) 

GpsmOc \(ms + mujuifj, J 

The T- violating muon polarization (Pt) comes from the interference of IB and 
the imaginary part of the SD. Only the pseudo-scalar (Gp) and right-handed current 
(Gp — Gv + Ga) terms contribute: 

Prix^y) = P^{x,v)+P^{x,y) (9) 
P^{x,y) = av{x,y)[Im{AA + Av)], 



with 



Pri^^y) = [(^v{x,y) - crA{x,y)]Im,{A 



P 



(10) 



where crv(x, y) and aA{x, y) are analytic functions of the Dalitz plot variables {x = 
2Ej/Mk, y = 2E^/Mk)- Figure |^ shows the contours of muon polarization along 
and perpendicular to the muon momentum within the decay plane, and the contours 
of cry and cry — cr^ as a function of the Dalitz plot variables, E^ and Ej. Events 
with high muon momentum have higher sensitivity to the T- violating parameters. 

The detectors designed for K'^ — > ^'^'k^v muon-polarization measurement can 
be used for studying /i+i^7 decays, with only minimal modification. In the 

near future, the first measurement of the T- violating polarization in ii'^vj is 

expected from KEK-E246. It is expected to have sensitivity of the order of 5% 
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Figure 1: Muon polarization for — > decays along the muon momentum 

(top left), normal to the muon momentum in the decay plane (top right), and ay 
(bottom left) and ay — oa (bottom right) on the Dalitz plot. The out-of-plane 
polarization would be proportional to the numbers on the contours in the bottom two 
plots, as given in Eq. [7^ 
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A — > event is usually identified in a detector through a coincidence of 

a stopped muon with an energetic acoUinear photon. To reject background from 
decays with tt*^ (mainly fi^n'^v), events with additional photons are 

rejected. This requires a hermetic photon veto system surrounding the decay 
region. With a precise measurement of the photon energy and the muon momentum, 
the K'^ — !■ decay can be fully reconstructed. This helps to further reduce 

background. A precision calorimeter with good spatial resolution is essential. The 
muon momentum can be measured in a magnetic field, as in KEK-E246, or the 
muon energy can be derived from the muon range in the polarimeter, as proposed 
in the new BNL experiment. A study from the latter showed that with sufficient 
shielding around the decay volume, and photon veto down to about 20MeV, one 
could obtain a sample of 3 x 10^ K'^ /i+i/'-v events, with a signal to background 
ratio of about 2:1, in 2000 hours of running.cEl This would provide constraints on 
the new interactions at sensitivities of: 

5lm{A^) = 20X10^3 

An experiment with the sensitivity of Eq. |ll|, will clearly probe new physics. Events 
with high muon energy have even higher sensitivity to T-violation; they are also 
beyond the kinematic limit of /i+Tr'^t^ events. They have a much better signal 

to background ratio, but at the expense of efficiency for signal. This is partially 
compensated by the higher sensitivity. For example, selecting events with high 
muon energy {2E^/Mk > 0.95), the average cry increases from 0.11 (for all events) 
to 0.15; the effect is larger for (cry — a a), which increases from 0.04 to 0.25. We 
note that this is also the region where the asymmetry from FSI may be maximal. 
The selection of events on the Dalitz plot can be optimized differently for a search 
for pseudo-scalar or right-handed currents. 

It should be noted that the T-conserving components of the muon polarization 
in K'^ — > ^+1^7 and —^ ^~^Tr^v decays are of opposite sign. By measuring the 
T-conserving polarization of the accepted —^ ^^v^ sample, the background 
from + fjL^n^v can be easily evaluated. 

4. Kl ^ ^Ji+^i- 

This decay has been measured recently with high statistics.S Experiment E871 at 
BNL has collected more than 6200 events with about 1 % background contamina- 
tion. The new branching ratio is: 

T{Kl TT+TT ) 

Using the measured branching ratio B{Kl tt+tt^) = (2.056 ± 0.033) x lO^'^ from 
Ref. H, we obtain B{Kl Ai+M") = (7.14 ±0.17) x m-^. The phenomenology of 
this reaction has been studied extensively in the past. I 
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The longitudinal muon polarization in this decay violates CP invariance. In 
general, the correlation between the /i+ and fi^ polarizations also contains infor- 
mation about CP violation and new physics, however this is much more difficult 
to measured The decay amplitude for Kl — > /^^M^ is known to be dominated by 
the two photon intermediate state. Interference of this amplitude with any other 
flavor-changing neutral scalar interaction can produce a longitudinal polarization. 
Following Herczegj£j we set the amplitude for Kl — > to: 

M{Kl ~> Ai^A^^) = a,u{p-)^^v{p+) + bu{p-)v{p+) 

a = a2 + i€ai (13) 
b = 62 + iibi 

where 02 and 62 are the CP conserving and CP violating amplitudes of K2, respec- 
tively. Similarly, ai and hi are the CP violating and CP conserving amplitudes for 
Ki, respectively. The decay rate and the longitudinal polarization (Pl) are then 
given by: 

r = !!^(H2 + ;32|^|2) 



= -^----H- ' (14) 

I„I2 1 q2\h\2 ^ ■> 



2/3Im(6a*) 
(3 = (1 - 4m2/TO|-)i/2 ^ 0.905 



The expression for the decay rate has several parts: terms of Oe^ can be ne- 
glected, terms of Oe multiply amplitudes with direct CP violation in Ki or K2 and 
are therefore small. The remaining terms are Re(a2)^ -I-Im(a2)^ -l-Re(62)^ -f Im(62)^- 
Each of these reflects a sum of contributions from electroweak interactions and pos- 
sible non-electroweak physics. The lowest-order standard-model electroweak con- 
tributions arise from the intermediate two photon diagram and the short-distance 
diagrams involving mainly the top quark, W, and Z exchanged 

The largest contribution to the total branching ratio is known to be the absorp- 
tive contribution (Im(a2)) from two photon exchange. @. This can be calculated in 
a model independent way: 



r(ifL^77) 'ruK' 2/3 V " 1 - 

Eq. |l^ is also known as the unitarity bound. All other contributions to Im(a2) are 
much smaller. Im(62) in the standard model is due to K2 tttt — > 77(CP=-|-1) —> 
this is constrained to be very small by e'/e. The short-distance electroweak 
diagrams contribute to Re(a2). The branching ratio from this contribution is also 
well calculated in the standard model to be: 



BsDiK^^^+p-) = ^^^i^^l^J^[Y,ReiX,) + Y,Re{\,)f (16) 

tk+ tt^ sm t'lvlKsl 
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Yq are Inami-Lim functions of Xq — M^/M^^; Xj — V*gVjd are combinations of the 
CKM matrix elemcnts.Ell Or, in terms of the Wolfenstein parameters p and A: 

Bsd{Kl ^ /i+M") = 1-51 X 10-^A4(po - pf (17) 

where p — p{\ — A^/2), po = 1-2 is the charm contribution. 

The second contribution to Re(a2) comes from — > 7*7* M^M i in which 
the photons are off mass-sheh. The sign and magnitude of this contribution are stiU 
uncertain; this uncertainty hmits the precision with which the experimental mea- 
surement can be used to obtain a hmit on polarization, as well as on the fundamental 
standard-model parameter p. First order amplitudes involving new flavor-changing 
vector bosons, leptoquarks, or scalar particles, such as a light Higgs, contribute to 
Re(62); these contributions would be responsible for a possibly large longitudinal 
polarization. 

The unitarity bound from the absorptive amplitude can be subtracted from the 
measured branching ratio to obtain the sum of all the remaining parts. In Ref. 
this subtraction is performed using the measured ratio T{Kl p^ p^)/T{Kl 
TT+TT^) and the unitarity bound calculated for the ratio T{Kl — > 'yj)/r{KL 
TT+TT^), with proper account of correlations in errorsEa They obtain a dispersive 
contribution to the branching ratio of (0.11 ± 0.18) x 10~^: 

^^((Re(a|^) ± Re(a^^))2 + /9'Re(52)2) = (0.11 ± 0.18) x 10"=' x T{Kl) (18) 
87r 

Here Re(a2^) is the magnitude of the amplitude due to the electroweak short- 
distance diagrams and Re(a2^) is the magnitude of the dispersive contribution 
from the two photon diagram. We will now consider the implications of this result 
on the value of the polarization: 

Concentrating only on the direct CP violating contribution, the longitudinal 
polarization can be written: 

= ~ A viK''^\ -, Re(j^2)Im(a2) (19) 

One can limit the magnitude of this polarization using the measurement in Eq. |l^, 
by assuming Re(a2) = 0. Then, for the square of the polarization, we can write: 

IP.P < 4 X (l - (20) 
V T{Kl p+p ) / T{Kl p+p ) 

Using the Particle Data GroupEl prescription, we obtain an upper limit" for the 
polarization of: 

\Pl\ < 0.44 90%C.L. (21) 

"Following the PDG prescriptioripthe value of IP^P is 0.058 ± 0.094. However, it should be noted 
that a previous result from E791m has an unphysical central value of —0.06 ± 0.21 for \Pl\'^- 
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Clearly, a better limit can be obtained if Re(a2) can be specified more accurately. 

There are several mechanisms that could cn n t r jb uto towards generating a signif- 
icant value of Re(62) and a large polarization.ll3li3-Lj'L3 In general, Pl is quite sensi- 
tive to the presence of light scalars with CP violating Yukawa couplings. The polar- 
ization can originate in the standard CKM model because of the induced s ~ d — H 
coupling, where H is the standard-model Higgs. Very large polarizations can be 
expected for a light Higgs boson with a mass comparable to Kl- Results from LEP 
that limit the Higgs boson massesto be larger than 115 GeV, rule out detectable 



polarization from this mechanism.cj In Ref. 63 multi-Higgs, leptoquark, as well as 
left-right symmetric models are considered. Polarizations of the order of few per- 
cent are possible, without violating bounds from neutron or electron electric dipole 
moments, or from lepton-flavor violation. Supersymmetric contributions to are 
expected to be small (~ 10~^), without fine-tuning of parameters. 

The dominant contribution to Pl within the standard model is through the 
indirect CP violation parameterized by e. The direct CP violating amplitudes (62 
and ai) can be neglected, and using the value of e:i@ 

e = 2.3 X lO^^e*'* (22) 
= 43.7° ± 0.1° (23) 

The polarization can be written as 

= --^TTIT-^^^^^^^— 7skl[Re(a2)(Re(6i)-Im(5i)) 
47rl [Kl H'^H )v 2 

+ Im(a2)(Rc(6i) + Im(fei))] (24) 

Ecker and Pich have calculated the value of bi using chiral perturbation theory 
(CHPT), to obtain a rather good prediction for PlI^ They obtain: 



(1.9)1.5 < \Pl\ ■ W'J-^^^—^^ < 2.5(2.6) (25) 

The numbers without (with) brackets correspond to 1 standard deviation, a, (2 
(t) errors. There are three sources of uncertainty in the above estimate. The un- 
certainty on the octet coupling strength (Gs) in CHPT is largely eliminated by 
normalizing to the measured B{Ks 77) = (2.4 ± 0.9) x 10^^.Ea The uncertainty 
in the upper limit is small because it comes from Im(a2). The uncertainty in the 
lower limit comes from the unknown sign and magnitude of Re(o2). Ecker and Pich 
point out that there is a constructive interference in the term multiplying Im(a2), 
which makes their estimate larger than a previous estimate.il To summarize, the 
polarization due to indirect CP violation is estimated to be ~ 2 x 10"'^ with 30-40% 
uncertainty. 

The main experimental difficulty in this measurement is the small branching 
fraction of (7.14 ± 0.17) x 10^^ for the decay. Much effort must therefore be put 
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into separating these events from background, before a polarization analysis can 
be performed. Experiment E871,E3 with ~ 6200 events, was optimized to look 
for Kl /i^e^. We have estimated that, if the experiment were optimized for 
Kl IJ-^f^ , and the beam intensity were increased, E871 could collect about 
20,000 events in two years of running. Appropriate upgrades to the marble muon- 
range detector will allow approximately 50% of the positive muon decays to be 
analyzed. A study by the E791 collaboration (previous version of E871), indicated 
that, with 20,000 events, they could achieve a sensitivity of 11% on the longitudinal 
polarization asymmetry of the positive muons.BEl 

The ~ 1% background to Kl — *■ /^"""/i" arises mainly from Kl ^ tt^ fjL^v events, 
in which the charged pion decays in flight or is misidentified as a muon, and the 
momentum of one of the charged particles is mismeasured so that the jiji invariant 
mass is higher than the kinematic limit of 489 MeV. The positive muons in the 
'K~ pL^u background events at the kinematic endpoint will be almost completely 
longitudinally polarized; the positive muons in the background — » j-i^v will 
come from tt"*" decay, and will be polarized (with strength dependent on the 7r+ decay 
angle within the experimental acceptance) in the opposite direction. Given these 
considerations, this background will not be problematic for an experiment with a 
sensitivity of ^ 0.1, however, a measurement of the — > M^/^^ polarization better 
than 1% will require less background. 

5. —f TT^ jl'^ fl^ 

This decay has a very rich structure that can lead to important measurements. 
Table ^ shows three different asymmetries that could be interesting to measure: 
longitudinal muon polarization, transverse muon polarization, and transverse fi^ 
polarization in correlation with longitudinal polarization. The decay has re- 
cently been analyzed extensively.E3"L3 The different processes that govern the de- 
cay are: the one-photon intermediate state, a two-photon intermediate state, the 
short-distance "Z-penguin" and "W-box" graphs, and potential contributions from 
extensions to the Higgs sector. 

The dominant amplitude from the one photon intermediate state (with a vector 
form factor) is best understood in the framework of chiral perturbation theory, 
although experimental inputs are needed to fully describe the decay.ES The 
TT^l^l^ decays can be discussed using the following variables: 



(Pi +P2) 



2 



m 



K 



y = 



2P-{p,-p2)_ 

) 

(26) 



where pi and p2 are the four-momenta of the negative and positive lepton. P 
and k are the four-momenta of the kaon and the pion, respectively. A(a, 6, c) = 
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+ 6^ + — 2(ab + bc + ac) and = m\/m\. The limits of the phase space are 
given by: 

4to^ ^ _ ^ {rriK — m-jr) 



< X < 



\y\ < (27) 



xm 



K 



In the kaon rest frame, the energy of the pion is given by: E^^ = (Afif /2)(1 — x — x,r)- 
In terms of these variables, the decay rate is given by 

""^ :\''\l,x,x.)\M\^ (28) 



dxdy 64(27r)3 

where M. is the matri -^ elem ent. For a vector interaction model, the decay rate is 



given by the following!^ 

dxdy IGtt 



X'/\l,x,x^){l-y^mx)? (29) 



where G% is the octet coupling constant. The general features of if ^ ■nl'^l~ decay 
spectrum can be understood in terms of angular momentum conservation. In the 
vector interaction model, the quantum numbers of the I'^l^ pair are J^*^ = 1 . 
The lepton pair must therefore be in a p state relative to the pion, and it must be 
longitudinally polarized. This has two consequences: the angle between the lepton 
and the pion in the rest frame of the lepton pair has a distribution that goes as 
(1 — cos^ 6*) (This is the reason for the (1 — y^) term; y m cos9, is the difference 
in the energies of the lepton pair divided by the pion momentum in the kaon rest 
frame); as the mass of the lepton pair gets large, the orbital angular momentum 
barrier forces the decay spectrum to fall faster than phase space, and to vanish at 
the kinematic endpoint. 

The form factor <j){x) has been measured very well by experiment E865, with 
more than 10,000 events in the K'^ — > T:~^e~^e~ decay modeJiJ The form factor 
is linear in x, with a large slope parameter that has been measured with high 
precision. The dominance of the vector current in the decay has been demonstrated 
experimentally. The measurement puts significant constraints on models containing 
long-distance contributions to if ^ 7rZ+Z~ decays. The A'+ 7r+/i+/i~ decays 
have been measured in two different experiments, ^ 200 events in E787,0 and 
~ 400 events in E865.0 The average branching ratio is 7.6 x 10~*, but the two 
experiments disagree by 3.3a. 

S(iC+ ^ 7r+/i+^~,E787) = (5.0 ± 0.4(stai) ± 0.7(sysi) ± 0.6(i/i)) x 10"®; 
B{K+ ^^^+^I+^i-',E865) = (9.22 ± 0.60(siai) ± 0.49(syst)) x 10"^ (30) 

Using the more recent if + 7r+e+e" measurement of the form factor, the if + 
7r+^+/i" branching ratio is predicted to be (8.7 ± 0.4) x 10"®. The E787 result 
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used an old value for the form factor; the theoretical error in their measurement 
corresponds to this form factor. The disagreement between the two measurements 
is large even if we correct the E787 result using the new form factor. 

In general, the decay process can have contributions from scalar, vector, pseudo- 
scalar and axial- vector interactions, with corresponding form factors, i^s, -Fy, -Fp, 
and Fa (following the notation in Ref. 

M = Fsu{pi,s)v{pi,s) + Fpu{pi,s)j5v{pi,s) 

+Fvp'^u{pi, s)'-ff,v{pi,s) + FAPkHPh s)'-ff,'j5v{pi,s) (31) 

Here pk, P-k, Pi, and pi are the kaon, pion, lepton, and antilepton 4-momenta. Any 
interference between the terms with a complex phase difference leads to polarization 
effects. Within the standard model, the largest contribution is from the one-photon 
intermediate state to the vector form factor [Fv], which is expected to be almost 
real. The scalar form factor is expected to get a small contribution from only 
the two-photon intermediate state. .Fp and Fa get contributions from the short- 
distance "Z-penguin" and "W-box" diagrams, where the dominant term arises from 
t-quark exchange, both form factors are therefore proportional to VtsV^^, with a 
small contribution from the charm quark. 

The parity- violating longitudinal muon polarization (asymmetry (4) in Table |l|) 
has terms proportional to Re(FpFy) and to Re(FyF'^). It is therefore sensitive 
to the Wolfenstein parameter p. The value of this polarization within the stan- 
dard model is estimated beyond the leading logarithms to be 0.003 — 0.0096, for 
— 0.25 < p < 0.25, and depends on the experimentally accessible region of phase 
space.cj Neglecting the dependence of the form factor, the polarization asymmetry 
varies as (1 — 4TO^/xm^)A"^^^(l, XjXt^)- It has a maximum at a /i/i invariant mass of 
approximately 250 MeV/(? . There is a non- negligible contribution to this polariza- 
tion from the long-distance two photon process, which has not as yet been calculated 
accurately.tj The T- violating transverse polarization (asymmetry (5) in Table |l|) is 
proportional to lTa{FsFy), and it is therefore expected to be small within the stan- 
dard model, and the final-state interaction correction to this polarization is expected 
to be ^ 10^'^. The T-violating spin-spin correlation that involves both /i+ and /i^ 
polarizations (asymmetry (6) in Table |l|) has terms proportional to ~lu\{FyFp) 
and Im(FvF^). This is expected to have much smaller final-state interaction cor- 
rections, and is theoretically sound. Within the standard model, such asymmetry 
is proportional to the CKM parameter -q, and it is expected to increase with 
invariant mass, becoming as large as ~ 0.06 in some parts of the decay phase space 
(see Fig. A good measurement would certainly be valuable in understanding 
CP violation. Similar to Kl /i^/i^, this asymmetry could also receive contribu- 
tions from non-standard physics, such as from charged Higgs or leptoquarks.E3 The 
consensus in the literature on the numerical results is summarized in Table |^. It 
should be noted that these estimates were made before the precise measurement of 



the form factor Fy in Ref. |70|, nevertheless, that should not affect substantially the 
values in Table |2[ 
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Figure 2: T-violating spin-spin correlation, or asymmetry (6) in Table mr-f or 
7r+/i+/i^ over the Dalitz plot. (After Agrawal, Ng, Belanger, and Geri^3). 



The main experimental difficulty is in isolating the rare —>■ 7r+/i+/i^ decays 
from background. For both E787 and E865, the main background is the 
7r+7r+7r~ decay, which has a branching ratio 6 orders of magnitude higher. This 
background must be suppressed at the trigger level as well as in the final analysis. 
In E865, the charged pions can be misidentified as muons due to pion decays in 
flight. In E787, where the unique total kinetic energy is used as the signature, the 
additional energy release in 7r~ nuclear capture can increase the apparent kinetic 
energy of tt+tt+tt^ events. Neither experiment was optimized for — s- 

TT^/i+zi", and both experiments took dedicated triggers in only a 

limited running period. 

Scaled from E865, a dedicated and optimized experiment, with improved ac- 
ceptance and trigger capability for muons, increased kaon flux with reduced beam 
halo, longer running time, and an excellent muon polarimeter, could yield 50, 000 
events. This could result in a measurement of /i^ polarization with 5Pl ~ 0.05, as- 
suming the same analyzing power as the estimate for the polarization measurement 
of Kl ^'^ . An interesting option for this measurement could be to rely on 
the decay-at-rest technique used by E787: the A'+ tt+tt+tt" background could 
be suppressed by selecting the tt"*" momentum to be above the kinematic limit of 
— > 7r+7r+7r~ decays (125 MeV/c); then both /i+ and jjT would have low kinetic 
energies, and could be stopped and analyzed in the same active target that serves to 
stop incoming kaons. Clearly, measuring asymmetries that require analyzing both 
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and /i^ polarizations will be difficult because ^~ decays have a lower analyzing 
power due to muon capture into atomic orbits around nuclei in the muon stopper. 

Authors of Ref. ^ point out that there are only a few interesting kaon mea- 
surements that offer clean theoretical interpretation in terms of the standard weak- 
interaction parameters. The measurement of the branching ratio of tt'^vv 
provides a rneasiirement of |VJ*Vi(i|, the branching ratio —^ t:^vv is propor- 
tional to 77^.l2rL3 The measurement of the polarization asymmetries (4) and (6) 
from Table |l| could be of comparable significance in terms of the overall theoreti- 
cal uncertainties. The experimental difficulty should be similar to other proposed 
investigations of the K and B systems. 

6. Kl ^ -K^fi+fi- 

The structure of 7r°?+Z~ decays is more complex than that of tt+Z+Z" 

decays because of CP suppression. There are three possible contributions to the 
decay amplitude: 1) direct CP-violating contribution from electroweak penguin 
and W-box diagrams, 2) indirect CP-violating amplitude from the Ki 7r°Z+Z^ 
component in K^, and 3) CP-conserving amplitude from the ir'^jj intermedialfi 
state. The sizes of the three contributions depend on the final-state lepton.HlIa 
The CP-conserving two-photon contribution to the electron mode is expected to be 
(1 — 4) X 10^^^, based on Kl 7r°77 data. Although suppressed in phase space, 
this contribution to the muon modejs comparable to the electron mode because of 
a term proportional to lepton mass.ES The CP-violating contribution to the electron 
mode is expected to be 



CPV 



15-3as - 6.8^as -1- 2.8 



X 10^12 (32) 



and the muon mode is suppressed by a factor of 5 due to phase space.E3 as is an 
unknown parameter in the Kg ~^ ir'^l^l^ vector form factor, and At = V^* Vt^. 

The modes 7r"e+e~ and Kl — > ir'^ fj,~^ fJ,~ have not as yet been observed; 

the current best limits on the branching ratios for Kl tt'^Z+Z^ were obtained by 
the KTeV experiment aLENAL; B{Kl n^'fi+fi-) < 3.8 x 10"^" and B{Kl 
TT°e~^e~) < 5.1 X 10~^°.cSlZI These limits were based on 2 observed events in each 
case, and expected backgrounds of 0.87 ± 0.15 for the muon mode and 1.06 ± 0.41 
for the electron mode. 

The main backgrounds for the muon mode were estimated to be from ^"'"/i~77 
(0.37±0.03) and 7r+7r-7r° (0.25±0.09), in which both charged pions decay in flight. 
Of these, the former background could be irreducible and therefore of great concern. 
The decay Kl /i^/i^77 proceeds via the Dalitz decay Kl ^ /i^/i^7, with an 
internal bremsstrahlung photon. KTeV has detected 4 such events, with expected 
background of 0.16 ± 0.08 resulting in a branching ratio of B{Kl /i'^/z~77) = 
(10.4l:5;g±0.7) X 10"^, consistent with the expectation of (9.1±0.8) x 10"^ obtained 
from QED and the measurement of Kl fJ-'^MT^i^ Such decays were analyzed 
by Greenlee as background to Kl 7r°e+e^.L3 He showed that the background 
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from e+e^77 can be suppressed by removing events with photon-photon invariant 
masses near the tt^ mass, and by constraining the energies and angles of accepted 
photons. Unfortunately, the criteria that can be used for suppressing e'^e~77 are 
not effective for the fi'^fi^j^ channel because the invariant mass of the photons 
is restricted by the energy available in the latter decay. In KTeV, the tt" mass 
resolution was expected to be 2.4 MeV; they apphed a 2.5 a cut around the 7r° 
mass to obtain the estimated background of 0.37 ± 0.03 from /i+/i^77. For any 
future experiment wishing to measure the polarization of the muon, a good range 
measurement of the muons will be required, which will most likely suppress the 
other backgrounds found in KTeV (7r"'"7r~7r'^ with pion decay to muons, and tt^/i^i/ 
with an accidental ir^). However, it seems unlikely that the background due to 
ji^ ^^^^ can be lowered. The single-event sensitivity of the KTeV result is quoted 
as 7 X 10~^^; the signal to background ratio, assuming that only /x"'"/x~77 will 
contribute in a future experiment, will therefore be around 1/5, if the standard- 
model signal is taken as B{Kl — > n'^ fi'^ n^) ^ 5 x 10^^^. 

The interesting direct CP-violating component must be extracted from any sig- 
nal found for Kl tt^I^I^ in the presence of two formidable obstacles: 1) the theo- 
retical uncertainty on contamination from indirect CP-violating and CP-conserving 
contributions and 2) the experimental background from the A complete 

discussion of this situation is beyond the scope of this article. In other reviews, it 
is stated that, in the case of —> n'^e'^e^ other measurements may be needed 
to understand the direct CP-violating contribution: e.g., the branching ratio for 
Ks — * TT^e^e^ and the energy asymmetry between e"*" and e^.L3c3 The branching 
ratio measurement of the tt/^/i mode provides additional information, such as an al- 
ternative measurement of the contribution due to the single-photon and two-photon 
intermediate states. The measurement of muon polarization in the irfifi case could 
also provide important additional information. This can be seen as follows: the 
decay amplitude can be divided into scalar (S), pseudo-scalar (P), vector (V), and 
axial-vector (A) parts. Of these, the scalar piece, which comes from the two-photon 
intermediate state is CP-conserving and the indirect CP-violating contribution is 
mostly vector. The short-distance direct CP violating amplitude has both vector 
and axial-vector parts. The longitudinal muon polarization (asymmetry (7) in Ta- 
ble while not strictly CP-violating, can arise only from interference of P and 
A with S and V, and therefore it provides information on the direct CP-violating 
amplitudes. Unlike for K~^, in the case of Kl decay, all the amplitudes should be 
of the same order of magnitude, and polarization effects should therefore be very 
large (~ 1), unless there are strong cancelations. The interference of scalar and 
vector components can give rise to transverse muon polarization (asymmetry (8} 
in Table This was studied in the context of a CHPT calculation of 0(p^)nJ 
The transverse polarization is displayed in Fig. |^ as a function of z = rnj^^/mj^, 
for several values of the renormalization parameter Ws- In this calculation, Im(ws) 
quantifies the direct CP-violation in the vector part of the matrix element. At the 
time of that study, the parameter Re(ws) was thought to be constrained by the 
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measurement of — > 7r+e+e^. Later, it was argued that the constiaint is model 
dependent, and Re(i(;s) should be considered as unknown as as-EilH Although a 
new detailed calculation is needed, especially to incorporate the new knowledge of 
the CP-conserving contribution, it is clear that the transverse polarization can be 
large. 

These large asymmetries should be easy to measure with sufficient statistics. 
This is no longer out of question, if one considers that the proponents of the BNL 
experiment KOPIO expect to measure ~ 50 Kl n'^vD events.Ej Measuring the 
muon polarization asymmetries in —> 7r°/i+/i^, together with the branching 
ratio and the lepton energy asymmetry, could be a good way of defeating the in- 
trinsic background from CP-conserving and indirect CP-violating amplitudes and 
the experimental background from /i^"/i^77. 

7. Conclusion 

Muon polarization from kaon decays has a rich phenomenology. In the case of 
Kl /i^/i^, —f TT^ ji^ ji'' and Kl 7r'^/i+/i^, new measurements could lead 
to important constraints on the CKM parameters, in particular the Wolfenstein 
parameters p and 77. The experimental difficulties should be comparable to those 
facing the rare kaon decay measurements of —>■ i:'^vT> and Kl —> ir'^vv, which 
are considered the best modes for understanding short distance physics in the kaon 
sector (for a recent review, see Ref. In particular, the flux available for the 
new experiment E926 for Kl — > t^'^vv- could be sufficient for a measurement of 
polarization in Kl 7r'^//+/i^.E3 As shown in Tables ^ and ||, for many cases, 
limits on the muon polarization will probe new physics beyond the standard model. 
In particular, the polarization will be sensitive to the physics of a more complicated 
Higgs sector, or leptoquarks, that could give rise to CP or T violation from sources 
outside of the standard model. 

We have examined the measurement of the out-of-plane muon polarization in 
K'^ —> ^~^Tr^u and — )■ fi'^i/j decays. Such measurements will not be sensitive 
to sources of CP violation in the standard model. Nevertheless, the measurements 
can be performed with sensitivity approaching SP ~ 10""* for —> ii~^Tr^v, and 
SP ~ 10^"^ for /f+ jjL^v^. For h^tt^v decays, this is well beyond the 

current direct limit of(— 4.2±4.9)xl0^'^, and the indirect limit of ^ 10^'', available 
from other experimental constraints. Although the electric dipole moments of the 
neutron and electron are considered more favorably for T violation outside the 
standard model, they do not cover the entire spectrum of possibilities beyond the 
standard model. At the moment, the measurement of T- violating polarization in 
/i+TT*^!/ and — > fx^v^ decays is well justified and should be considered 
complementary to other efforts in understanding CP violation. 

We thank Laurence Littenberg, Larry Trueman, Yoshi Kuno, Steve Kettell for 
useful discussions. This work was supported by DOE grant DE-AC02-98CH10886. 
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Figure 3: Up-down asymmetry in Kl — » 7r''/i+/i^ decay as a function of z (defined 
as z = ni?f^^/m\), for two possible values of Re{ws) and three different values of 
Im{ws) covering the expected ranges for this parameter. The different curves cor- 
respond to the following values of (ReWs, ImWs): double- dot- dashed (0.73, —10~^), 
long-dashed (0.73, 0), dot-dashed (0.73, +10-^), dashed (-1.00, -IQ"^ ), solid (- 
1.00,0), dotted (-1.00,+lQ-^ ), (after Ecker, Pich, and de RafaelB) 
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